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ABSTRACT 

Using moderate-resolution optical spectra from 58 background Lyman-break galaxies and quasars 
at ^ ^ 2.3 — 3 within a 11.5' x 13.5' area of the COSMOS field 1200 deg“^ projected area density or 
^ 2Ah~^ Mpc mean transverse separation), we reconstruct a 3D tomographic map of the foreground 
Lyct forest absorption at 2.2 < z < 2.5 with an effective smoothing scale of lsd ~ 3.5 Mpc 
comoving. Comparing with 61 coeval galaxies with spectroscopic redshifts in the same volume, we 
find that the galaxy positions are clearly biased towards regions with enhanced IGM absorption in the 
tomographic map. We find an extended IGM overdensity with deep absorption troughs at 2 : = 2.45 
associated with a recently-discovered galaxy protocluster at the same redshift. Based on simulations 
matched to our data, we estimate the enclosed dark matter mass within this IGM overdensity to be 
= 2.45) = (9 ± 4) X 10^^ h ~^ Mq, and argue based on this mass and absorption strength that 
it will form at least one 2 : ^ 0 galaxy cluster with M{z = 0) = (3 ± 2) x 10^^ h~^ Mq, although 
its elongated nature suggests that it will likely collapse into two separate clusters. We also point 
out a compact overdensity of six MOSDEE galaxies at 2 : = 2.30 within a r ^ lh~^ Mpc radius 
and Az ^ 0.006, which does not appear to have a large associated IGM overdensity. These results 
demonstrate the potential of hya forest tomography on larger volumes to study galaxy properties 
as a function of environment, as well as revealing the large-scale IGM overdensities associated with 
protoclusters or other features of large-scale structure. 

Subject headings: cosmology: observations — galaxies: high-redshift — intergalactic medium — 
quasars: absorption lines — galaxies: clusters: general — techniques: spectro¬ 
scopic 


I. INTRODUGTION 

The study of high-redshift {z > 2) galaxy protoclus¬ 
ters is a topic of increasing interest, providing a route to 
studying galaxy evolution, AGN, plasma physics, and our 
models of gravity. They are also critical laboratories for 
understanding the growth of massive galaxy clusters at 
2 : ^ 0. While many studies have successfully found pro¬ 
toclusters around active supermassive bl ack holes such 


as radio galaxies or luminous quasars (see Hennawi et al 
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2QI5|and|GQQke et al.|[Ml4] fQr recent examples, or Ghi- 


ang et ahj 2013 for a compilation), it is difficult to tell 


whether these 'signpost’ protoclusters are representative 
of the overall population - which, at least in simulations, 
shows a wide range of properties. 

More uniform, ‘blind’, search techniques are therefore 
required for unbiased samples of protoclusters. Several 


tometric or galaxy redshift surveys (e.g., Steidel et al. 

2005 

Gobat et al.||201I| Gucciati et al.||2014l |Yuan et al. 

2014 

Casey et al.||20I5 

. Searching photometric redshift 

catalogs (jChiang et al. 

2014) is another promising new 


wavelength photometry that enable accurate photomet¬ 
ric redshfits. Unfortunately we expect protoclusters to 
be rare, and surveying large volumes of space with these 
techniques i s expensive. 

Recently, Stark et al. (20I5b|) (hereafter SI5) argued 
that 3D toniographic reconstructions of the intergala ctic 
medium (IGM) Lyman-a (Lyg) forest absorption (|Pi-| 
chon et al.||200It fCaucci et al.|[20081 |Lee et al.||2U14a|) 
can be used to efficiently search for proto clusters! A 
protocluster’s Lya absorption signature extends over a 
large region (r > Mpc) allowing it to be mapped 

by background galaxies separated by several transverse 
Mpc, corresponding to relatively bright limiting magni¬ 
tudes (< 24.5 mag) and therefore accessible to existing 
8-IOm class telescopes. In fact several studies have al¬ 
ready found Lya absorption associated with protoclus- 
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ters either in single sightlines (Hennawi et al. 2^^ or 
by st acking multiple background spectra ( [Cucciati et al 


2Q14[ Hayashino et al. in prep.). Moderate-resolution 


spectra with careful sightline selection should enable 3D 
map ping even with mod est signal-to-noise ratios (S/N). 


In Lee et al. (|2014b) (hereafter L14b), we made the 
first attempt at Lya forest tomography using a set of 24 
background Lyman-break galaxy (LEG) spectra within a 
5' X 12' area in the COSMOS field. This marked the first- 
ever use of LBGs, instead of quasars, for Lya forest anal¬ 
ysis, and the resulting reconstruction at 2.20 < 2 : < 2.45 
was the first 3D map of the z > 2 universe probing Mpc- 
scale structures. In this paper, as part of the pilot ob¬ 
servations for the COSMOS Lyman-Alpha Mapping And 
Tomography Observations (CLAMATO) survey to map 
the high-redshift IGM, we extend the L14b map volume 
by a factor of ^ 2.5 x. This volume starts to approach 
that where we expect to see protocluster candidates in a 
blind survey: in N-body simulations at ^ ^ 2.5 the co¬ 
moving number density of halo progenitors which hy z = 
0 have masses above 10^^ h~^MQ is ^ 10“^ Mpc“^, or 
to 10“^ Mpc~^ for the proge nitors of 3 x 10^^ h~^MQ 
halos (e.g. Stark et al. 2015b). In the simulations of 
S15, mock observations with similar sightline density and 
noise as CL AMATO found more than 3/4 of the proto¬ 
clusters which would eventually grow to form clusters 
more massive than 3 x 10^^ h~^MQ at z = 0 but only 
1/5 of the protoclusters which would form 10^^ 
halos hy z = 0. Our observations probe a volume of 
V = b.S X 10^ h“^Mpc^ at 2.2 < z < 2.5. To the extent 
that these simulations are accurate, there is a ^ 10% 
chance we should find a protocluster in the surveyed vol¬ 
ume (see also 

Two other developments in the past year lend addi¬ 
tional s ynergy: (i) the first data release of the MOSDEF 


survey (Kriek et al. 2015, hereafter K15), whose near- 
infrared (TORJTieEurax^ redshifts reduce systemic red- 
shift uncertainties hence allowing a cleaner comparison 
with our tomographic map; and (ii) the discovery of a 
2 : = 2.45 galaxy protocluster w ithin our target field (Di- 


ener et aL]|2015| hereafter D15; |Chiang et al.]|2015| here- 
after C15). 

As we shall see, our tomographic map suggests a com¬ 
plex 3D multi-pronged structure for this protocluster, 
while also revealing other intriguing structures. We also 
find that galaxies co-eval with our map preferentially in¬ 
habit higher absorption regions of the IGM, suggesting 
that Lya tomography does indeed probe the underlying 
large scale structure. 

This paper is organized as follows: we describe the data 
in Section [2 and tomographic reconstruction of the fore¬ 
ground IGM in Section and then describe simulations 
used for the subsequent analysis (Section]^. Section]^ 
then describes the protoclusters and overdensities founa 
within our map volume. 

In this paper, we assume a concordance flat ACDM 
cosmology, with Qm = 0.27, Oy = 0.73 and Hq = 
70 km s“^ Mpc“^. 

2. OBSERVATIONS 

We targeted 2 .3 < 2 : < 3 galaxies and AGN within the 
COSMOS field (Scoville et al. 2007 Capak et al.||2007) 
in order to probe the foreground IGM Lya forest absorp¬ 



Wavelength (angstroms) 

Figure 1. Examples of galaxy spectra used for our tomographic 
reconstruction (smoothed by a 3-pixel boxcar for clarity), along 
with the estimated pixel noise in red, with masked pixels set to 
zero. The Lya wavelength at the galaxy redshift is indicated by 
the vertical dashed-lines , while the green cur ves shows the LEG 
composite template from|Shapley et al.|(|2003|) for comparison. All 
these background sources probe the known 2 : = 2.44 protoclus¬ 
ter, with the pink shaded region highlighting the Lya absorption 
wavelength (A 4191 A) at the protocluster redshift. Strong as¬ 
sociated absorption is apparent even in these individual spectra. 
Orange curves indicate the estimated continuum level for each ob¬ 
ject. Note that object #00958 (3rd panel from top) has a broad 
Lya emission line characteristic of an AGN, but shows strong ISM 
absorption redwards of Lya and is thus treated like a galaxy for 
continuum-fitting. 


tion at 2 : ^ 2.3. Top priority was given to o bjects with 


confir med redshif ts from the zCOSMO S-Deep (Lilly et al. 


2007) and VUDS (Le Fevre et al. 2015) spectroscopic sur 
veys, while we also, where available, added grism redshift 
information kindly p rovided by the 3D-HST team (e.g. 


Brammer et al.| 

2012). We also selected targets based on 

photometric red 

[shifts from 

Ilbert et al. 

(2009) as well as 

Salvato et al. (2011 

) for X-ray detectec 

1 sources, which 


rate. This target selection process was c arried out in 
Febru ary 2( 314, prior to the publ ications by |Diener et al 


(2015) and Chiang et al. (2015) reporting galaxy over- 
densities in the field. 

We observe d these targets with the LRIS Dou ble- 
Spectrograph (|Oke et al. 1995 Steidel et al. 2004) on 
the Keck-I telescope at IVIaunakea, Hawaih, during 2014 
March 26-27/29-30 and 2015 April 18-20, using multi¬ 
object slitmasks. We used the B600/4000 grism on the 
blue arm and R400/6000 grating on the red with the 
d560 dichroicp^ although we work only with the blue 
spectra in this study. With our 1" slits, this yields 


R = A/AA ~ 1000, which corresponds to a line-of-sight 
FWHM of 3.2 h~^ Mpc at 2 : ^ 2.3. Despite inclement 
weather, we managed ^ 15hrs on-target, observing 7 


The first night of 2014 observations were taken with the 
R600/7500 grating and d500 dichroic. 
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Xperp(z=2.35) (h'^Mpc) 

15 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 -1 



Right Ascension 

Center: R.A. 10 00 17.95 Dec +02 15 50.4 


Figure 2. Hubble ACS image of the COSMOS field (|Koekemoer et al.|2007|) showing the position of our background galaxies and quasars 
on the sky, labeled by their ID number. The overall blue box represents the transverse extent of our tomographic map, while the seven 
narrow slices indicate the footprint of the map projections shown in Figure ^ The upper- and right-axis labels denote the transverse 
comoving distance at 2 : = 2.35, the mean redshift of the map. 


overlapping masks with 7200-9000s exposure times. The 
observations were mostly carried out in 0.5 — 1" seeing, al¬ 
though one mask was observed in « 1.4" seeing and con¬ 
sequently had a low yield of useable spectra. The data 
was reduced with the LowRedux package in the XIDL 
suite of data reduction softwarj^ and objects targeted 
in more than one mask had their spectra co-added. We 
extracted 162 unique spectra, of which 58 objects had 
secure visual identifications at redshifts 2.3 < 2 : < 3.0 
and therefore Lya forest absorption covering our desired 
2.2 < 2 ; < 2.5 range, as well as adequate S/N (S/N > 1.3 

http://www.ucolick.org/~xavier/LowRedux/lris_cook. 

html 


per pixel) within the Lya forest. The majority of these 
sources are LBGs identified through Lya emission or in¬ 
trinsic absorption, although our sample also includes 4 
QSOs, one of which is our brightest object {g = 20.12). 
The faintest object satisfying our S/N criterion was a 

= 25.18 LEG, although this was observed in two masks; 
the faintest object from a single mask has g = 24.85. 

Several example spectra are shown in Figure while 
the source positions are shown in Figure Tne lat¬ 
ter also indicates the 11.8' x 13.5' footprint of our to¬ 
mographic map. Within the map footprint, we have 52 
background sources which translates to a projected area 
density of ^ 1200 deg~^. An additional six sources lie 
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just outside the footprint but will be included in the re¬ 
construction input. The mean transverse separation be¬ 
tween the Lyof forest sightlines probing the foreground 
IGM is {di_) ^ 2.4 Mpc. 

We next estimate the intrinsic continua, C, of the 
sources. For the quasars, we apply PCA-based me an- 
flux regulation (MF-PCA; e.g., Lee et al.]2Q12 , 2013a) to 
the restframe 1041 A < A < 1185 A Ly g fore^ region us¬ 
ing templates from [Paris et al. (2011), masking intrinsic 
broad absorption where necessary. A similar process is 
applied on the gala xies, albeit assumin g a fixed contin¬ 
uum template from Berry et al.l (20121) and adopting a 


more generous Lya forest range (1040 A < A < 1195 A). 
We also mask ±5 A around possible intrinsic absorption 
at N II A1084.0, N I A1134.4, and C III A1175.7. L14b 
estimated an rms error of ^ 10% for the continua, which 
is adequate considering our noisy spectra. 

We then compute the Lya forest fluctuations, 6f, from 
the observed spectral flux density, /: 


6f 


f 


C{F{z)) 


1 , 


( 1 ) 


where (F(z)) is the mean Lya transmission from 
Faucher-Giguere et al. (2008) evaluated at redshift 
I'he dF values from our background sources, along with 
the associated pixel noise errors (Jn estimated by the re¬ 
duction pipeline, constitute a sparse sampling of the IGM 
Lya forest absorption, which we will reconstruct in the 
next section. 


3. TOMOGRAPHIC RECONSTRUCTION 

For the tomographic reconstruction, we define a 3D 
output grid with cells of comoving size 0.5h“^ Mpc, 
spanning a transverse comoving length at 2 ; = 2.35 of 
14 h~^ ]^c in the R.A. direction and 16 h~^ Mpc in Dec 
(Figure!^. The output grid spans 2.2 < < 2.5 along 

the line-of-sight, corresponding to 260 Mpc with our 
simplification that the comoving distance-redshift rela¬ 
tionship is evaluated at fixed 2 : = 2.35. The overall 
comoving volume is U = lAh~^ Mpc x 16Mpc x 
260 Mpc = 58240 h“^Mpc^ « (38.8 Mpc)^, i.e. 
nearly 2.5 x greater volume than the L14b map, which 
corresponds approximately to Slices #5-7 in Figure 
As in L14b, we use a Wiener Altering implementation 
developed by S15. The reconstructed flux field on the 
output grid is given by 

= Cmd • (Cdd + N)“^ • 8f, (2) 

where Cdd + N and Cmd are the data-data and map- 
data covariances, respectively. We assumed a diagonal 
form for the noise covariance matrix N = Na = •, and 

assumed a Gaussian covariance between any two points 
ri and r 2 , such that Cdd = Cmd = C(ri,r 2 ) and 


C(ri,r2) = cr|.exp 


(Ar.y 

2Ll 


exp 


2L\ 


(3) 


where Ary and Ar^ are the distance between ri and 
along, and transverse to the line-of-sight, respectively. 
Again following LI4b, we adopt transverse and line- 
of-sight correlation lengths of L±_ = 3.5 Mpc and 


I/|| = 2.7 h~^ Mpc, respectively, as well as a normaliza¬ 
tion oi (7F =0.8. These values were found to give reason¬ 
able reconstructions on simulated data sets with similar 
sightline sampling and S/N as our data. 

The resulting tomographic reconstruction of the fore¬ 
ground IGM is presented in Figure as a 3D visualiza¬ 
tion, while Figure shows the same map as a series of 
2D projections across Ax = 2 h~^ Mpc along the R.A. 
direction; a movie of the 3D visualization can be viewed 
onlinj^ As in L14b, we see large overdensities and un¬ 
derdensities spanning > 10 h~^ Mpc. The enlarged map 
volume also allows us to repeat the comparison of the 
map with coeval galaxies within the same volume. Using 
an internal compilation of spectroscopic redshifts within 
the GOSMOS held (albeit updated from the one used in 
L14b), we And 6 1 coeval galaxies, primarily from the 
zCOSMOS-Deep (ILilly et al.l [20071), VUDS (|Le Fevre 
MOSDEF (K15) sL 


et al. 


[2015 ) and M(JbDE4' (Ki5) surveys. We overplot 


these galaxy positions on Figure [^ These aforemen¬ 
tioned surveys are known to be uniformly flux-limited; 
we momentarily leave out galaxies targeted specifically 
as protocluster members, which will be discussed in Sec¬ 
tion O 

We compare the distribution of map values at the 
galaxy positions with the overall map PDF in Figure 
There is a clear skew in the absorption at the galax¬ 
ies’ positions towards negative 6^^ values, confirming 
that both the tomographic map and galaxy positions 
sample high-density regions. A two-sample Kolgomorov- 
Smirnov test finds P = 5.7 x 10“^ that the galaxy values 
are drawn from the overall distribution; this is a consid¬ 
erable improvement from the analogous comparison with 
18 coeval galaxies in LI4b, which had P = 0.24. 

There are several effects, however, that can introduce 
scatter into this comparison between galaxy positions 
and the large-scale IGM. Even in the absence of map 
reconstruction errors and with perfect redshift determi¬ 
nations, one expects a slight mismatch from the slightly 
different velocity fields traced by galaxies and the Lya 
forest, which is in turn due to the different halo biases 
of both tracers. However, in our case we expect the pri¬ 
mary sources of scatter to be reconstruction errors in 
the tomographic map, as well as line-of-sight positional 
uncertainties {ay ^ 300 km s“^ or aios ^ 3.5 h~^ Mpc) 
on the galaxies comparable to our map resolution. The 
latter arises from redshift errors induced by low spec¬ 
tral resolut ion (P < 200) ([Lilly et al. 2007 Le Eevre 
et al.[[2013) as well as intrinsic scatter between the true 
systemic redshift compared with redshift esti mates from 
Lya emission or UV absorption lines (e.g. [Ade lberger 
et al.[[200'5 Steidel et al.[[2010 Rakic et ^ 201l[ ). L14b 
showed using simulated reconstructions that this redshift 
error, along with reconstruction noise in the tomographic 
maps, would scatter some galaxies from overdensities 
into adjacent underdensities in the tomographic map. 

However, galaxy redshifts measured in the NI R have 
much smafl er redshift errors {ay r\j 60 km s [Steidel 


et al.|201Q|) due bo th to the higher resolution {R A 350 O 


Kriek et al. 2015) and the fact that rest-frame optical 
nebular lines are much better tracers of the systemic 
frame, so when we restrict the PDE comparison only to 


https://youtu.be/KeWlUJOPMYI 
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Figure 3. Three-dimensional visualization of the tomographic reconstruction from our Lyo; forest data, which can be viewed as a video 
online (24MB, https://youtu.be/KeWlUJOPMYI). The x- and y-dimensions, spanning 14 Mpc and 16 h~^ Mpc respectively, are on the 
sky plane (c.f. Figure!^, while the z-axis spans 260Mpc along the line-of-sight. All axes are in units of h~^ Mpc comoving, with 
the 2 ;-axis originating ^ z = 2.20 and extending to 2 : = 2.50. Since this map traces Lya forest transmitted fraction, more negative values 
correspond to higher-density regions. In the video, first 5 seconds show only the sightline positions before the actual map is revealed. 


Op ~ 

of 


the 31 galaxies observed in the NIR by MOSDEF, we find 
an even stronger skew towards higher map over densities: 
the MOSDEF galaxies sample a median map value of 
—0.24, whereas the other galaxies have a median 
~ —0.09. This comparison is encouraging and 
indicates that with expanded map volumes, we will be 
able to study correspondingly larger samples of coeval 
galaxies as a function of large-scale environment, with 
cuts in various galaxy properties such as star-formation 
rates, metallicities, AGN activity etc. 

It is worth briefly mentioning the large regions in our 
tornographic map with positive (blue regions in Fig¬ 
ure Q that are completely empty of galaxies. These likely 
co rrespond to the hig h-redshift voids recently discussed 


m 


Stark et al. (2015a), but we defer a more detailed in¬ 


vestigation to a subsequent paper. 


4. SIMULATIONS 

Before we proceed to search for overdensities in the 
COSMOS tomographic map, we first describe the N- 
body simulations that we will use to facillitate our sub¬ 
sequent analysis, which are the sa me simulation s as S15. 
These are 2560^-particle TreePM (White 2002) simula¬ 


tions with 8.6 X lO^h“^M0 equal-mass particles within 
a (256 Mpc)^ comoving volume. The assumed cos¬ 
molog y was consistent with |Planck Collaboration et al.| 
(|2013fc: ~ 0.31, 0.022, h = 0.6777, n, = 

0.9611, and erg = 0.83. 

Using a friends-of-friends algorithm to identify halos, 
S15 found 425 halos at 2 : = 0 with M > 10^^ h~^ Mpc 
that were defined as gaiaxy ciusters. The progenitor ha- 
ios for each ciuster were traced backwards to ^ = 2.5, and 
their center-of-mass at that epoch defined as the corre¬ 
sponding z = 2.5 protociuster positions. 

The veiocities and particies at 2 : = 2.5 were used to 
generate Lya forest spectra with the fluctuating Gunn- 


Peter son approximation (Croft et ai. 1998 Rorai et ai. 


2013). From these skewers, 815 also generated tomo- 


graphic maps from realistic mock data sets with vari¬ 
ous sight line sampling and noise properties. These re¬ 
constructions used the same Wiener-Altering code de¬ 
scribed in Section In this paper, we will work pri¬ 
marily with their simulated reconstructions with average 
sightline separation (d^) = 2.5 h~^ Mpc binned on to 
a {lh~^ Mpe)^ grid, which were designed specifically to 
match the L14b data set, including the spectral resolu¬ 
tion and S/N distribution in the mock spectra. 

Since progenitors of massive 2 : ^ 0 clusters oc cupy 
overdensities on ^ 5 — 10 h~^ Mpc scales at 2 : > 2 (Chi 


ang et al. 2013), S15 argued that an efficient way to 
identify ^ ^ 2.5 protoclusters in Lya forest tomographic 
maps is to first smooth the map with a cr = 4h“^ Mpc 
Gaussian kernel to obtain the smoothed flux and 
then identify regions with < —3.5 asm, where 
is the map variance in the respective map after smooth¬ 
ing. We And = 3.5 x 10“^ in our real map, and 
= 3.8 X 10 ^ the simulated maps. 

In Figure we show the distribution of smoothed 
values from the overall (256 Mpe)^ simula¬ 
tion volume (dotted-blue histogram), while the red his¬ 
togram shows the minimum d|U/crsm within Ah~^ Mpc 
of each of the 425 protoclusters within the simulation; 
the shaded histogram shows the subset of 27 protoclus¬ 
ters that will grow into M > 5x10^^ h~^ Mpc clusters by 
2 : = 0. This clearly shows that protoclusters occupy ex¬ 
treme Sp^/asrn regions that are well-separated from the 
overall distributior£3 


S15 showed a similar comparison in their Figure 3, but with 
the True’ Lya absorption field from the simulation, whereas we use 
the simulated tomographic reconstruction with noise and resolution 
matched to real data. 
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Figure 4. Projections of our Lyo: forest tomographic map across 2h~^ Mpc slices in the R.A. direction, with the slice locations indicated 
in Figure 1^ Note that negative values of 5^^^ (red colors) correspond to higher overdensities, while the white horizontal lines represent the 
skewer sam pling from the ba ckground sources. Grey squares mark the positions of coeval galaxies with spectroscopic redshifts from the 
zCOSMOS (|Lilly et al.|2007|) and MOSDEF (K15) surveys. The error bars indicate the LOS uncertainties associated with LEG redshift 
determination, while the galaxies without error bars have NIR redshifts, in which case the symbol widths denote the LOS uncertainty. 
Black contours in Slices #3-5 (at 2 ; 2.435 — 2.450) delineate regions with < —3.5f7sm, which is the protocluster criterion defined by 

S15. Pink circles in Slices #3-4 indicate the 2 : = 2.300 compact overdensity of MOSDEF galaxies within our volume. 


S15 have shown that this method of selecting proto¬ 
clusters yields good purity: applying the same selec¬ 
tion criteria on the simulated tomographic map with 
{d_\_) = 2.5 Mpc, we found 89 candidate protoclus¬ 
ters of which 81 will eventually collapse into a M > 
h~^ h~^ Mq halo by ^ = 0 90% purity); an¬ 

other 7 candidates 8%) have group-sized descendants 
of 3 X 10^^ Mq < M{z = 0) < 10^^ Mq, while 
only 1 candidate (^1%) will evolve into a low-mass halo 
{M{z = 0) < 3 X 10^^ h~^ Mq). These results are quali¬ 
tatively similar to those found by S15 but are not numer¬ 


ically identical despite being applied to the same simu¬ 
lations; this is likely because we have slightly different 
linking criteria for identifying contiguous thresholded re¬ 
gions. 

Figure shows several of the simulated protoclus¬ 
ters, in each case comparing three tomographic recon¬ 
structions matched to the L14b data, but with differ¬ 
ent random realizations of the sightline positions and 
pixel noise. The progenitors of massive clusters occupy 
such extended overdensities that they can be robustly 
detected with our sightline sampling, although at lower 
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Figure 5. PDF of the reconstructed absorption from our map 
(black), along with the corresponding PDF at the position of the 
61 coeval galaxies (red); negative absorption corresponds to higher 
overdensities. The galaxies PDF is clearly skewed towards over- 
dense regions of the map; this is even clearer in the case of the 
subsample of 31 MOSDEF redshifts (shaded red histogram). 


Data Reconstruction 
Sim Reconstruction 
Sim Protociusters 



Figure 6. PDF of the quantity S^^/asm used to select for pro¬ 
tociusters. The black histogram shows the distribution from our 
real map, while the blue-dotted histogram shows that from simu¬ 
lated reconstructions. The distribution corresponding to the 425 
M{z = 0) > 10^"^ Mq protociusters within the simulation is shown 
by the red histogram, while the hashed red histogram shows the 
subset of 27 protociusters with M( 2 ; = 0) > 5 x 10^^ Mq. The 
dotted-dashed red histogram shows the protoc luste r distribution 
including our toy feedback model (see Section [5.2| ). The shaded 
grey region indicates the S^^/asm < —3.5 criterion for selecting 
protociusters. Downward arrows indicate the values correspond¬ 
ing to the 2 ; = 2.435 and 2 : = 2.450 lobes of protocluster candidate 
found in our tomographic map, as well as the 2 ; = 2.300 overdensity 
of MOSDEF galaxies. 

masses {M{z = 0) ^ 10^^ h~^ M©) the absorption sig¬ 
nature can be less pronounced and is less likely to be 
selected by our criterion (e.g. bottom panel of Figure 
The overall completeness of the S15 selection criterion 
is ^ 20% for all M{z = 0) > 10^^ cluster progenitors 
(81/425), but this improves dramatically for higher-mass 
progenitors: we detect 41/51 protociusters with M{z = 
0) > 3 X 10^^ Mq 80% completeness) and 24/27 


protociusters with M{z = 0) > 5 x 10^^ h ^ Mq 90% 
completeness). The comoving space density of M{z = 
0) > 10^^ h~^ Mq clusters is n ^ 2.5 x 10“^ Mpc“^ 
within our simulation, which translates to an expectation 
of ^ 1.5 protociusters within our present observed map 
volume {V ^ 5.8 x 10^ h~^ Mpc). 

Figure also shows the smoothed S^^/asm distribu¬ 
tion from our COSMOS tomographic map (black his¬ 
togram), with both the dsm values derived individually 
for the data and simulated maps. There is some disagree¬ 
ment in the distributions from the data and the simula¬ 
tions, but the Lya forest absorption PDF is known to 
be a notoriously difficult quantity to model: our simu¬ 
lation used a DM-only N-body code that models only 
gravitatio nal clustering, but n ot the hydrodynamics of 
the ICM (White et al. 2010), nor have we attempted 
to accurately model systematics such as the continu um- 


errors, spec tral noise, optically-thick absorbers (e.g., Lee 


et al. 2015) , temperature-den sity relationship, or Jeans 
smoothing ( Rorai et al.||2013 ). For the present analysis, 
it is of sonie concern that our map PDF is somewhat 
larger at < —3.5 asm than the simulated PDF, but 
this is likely to be due to damped Lya absorbers (DLAs) 
and Lyman-limit systems (LLS’s) with column densities 


of Nm > 10^' cm“^ contaminating individual sightlines. 
By adding individual fake DLAs to our simulated re¬ 
constructions, we found that these could cause small 
regions with limited transverse extent (< lh~^ Mpc) 
to cross our protocluster threshold, but they cannot 
cause extended transverse IGM signatures characteris¬ 
tic of massive protociusters (Figure [^, since this re¬ 
quires strong correlated absorption across multiple ad¬ 
jacent sightlines. Over the total Az ^12 redshift path- 
length of our Lya forest sp ectra we expect only ^ 2 DLAs 


with A^hi > 10^^ '^ cm“^ (Prochaska et al. 2005), which 
makes it very unlikely that chance alignments of DLAs 
could cause a spurious protocluster candidate. Careful 
forward-mo deling of such systematics would be desirable 
to characterize marginal protocluster candidates in fu¬ 
ture data sets, but as we shall see we have only one pro¬ 
tocluster candidate in our data, which is highly extended 
(> 10 h~^ Mpc) in both the transverse and line-of-sight 
dimensions. 

S15 showed that various quantities related to the Lya 
forest absorption signature of simulated z = 2.b proto¬ 
ciusters can be related to their z = f) descendant masses, 
although they showed this mostly as a function of the 
‘true’ Lya absorption field in the simulation. We now 
carry out a similar analysis, but incorporating forward 
modeling by using the 81 protociusters identified through 
the tomographic maps of the same simulated absorption 
field, which incorporate realistic random sparse sampling 
of the intervening sightlines as well as noise and spectral 
resolution matched to our data. This can then be di¬ 
rectly co mpa red with overdensities identified in our data 
(Section 5.1). In the following discussion, ‘protoclus¬ 
ter candidates’ refer to 5^^ < — 3.5asm regions in the 
a = 4h“^ Mpc Gaussian-smoothed tomographic maps. 
Individual protocluster candidates are defined through a 
Mpc linking length regardless of physical continu¬ 
ity, so it is possible for a single candidate to be comprised 
of multiple disconnected structures, so long as their vox¬ 
els can be linked to within Mpc. 
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Figure 7. Simulated 2 ; = 2.5 galaxy protoclusters showu as projectious ou the sky plaue with thickuess 5h ^ Mpc aloug the liue-of- 
sight. Each map is ceutered ou the ceuter-of-mass of progeuitor halos that will eveutually collapse iuto 2 ; = 0 clusters with (top to bottom) 
M = [1.1 X 10^^, 8.8 X 10^"^, 1.2 X lO^'^] h~^ Mq. The left pauel ou each row shows the uuderlyiug DM distributiou (smoothed ou 1 h~^ Mpc 
scales), while the other three pauels show tomographic recoustructious from differeut realizatious of mock Lyo: forest data with similar 
sightliue separatious ({d±) = 2.5 h~^ Mpc) aud uoise as our data. The black coutours overplotted ou the tomographic maps iudicate 
regious that satisfy our protocluster criteriou of < —3.5crsm after smoothiug. While this threshold robustly selects massive cluster 
progeuitors, it is margiual for lower-mass protoclusters such as the bottom pauel. Note that iu the middle pauel, the lobe to the upper left 
([Ax, Ay] ^ [—10,8] h~^ Mpc) will iu fact a collapse iuto a separate, albeit lower mass (M = 1.5 x 10^"^ h~^ Mq) 2 ; ~ 0 cluster thau the 
ceutral object. 


Since the Lya forest absorption is a tracer of the un¬ 
derlying dark matter (DM) density, it should be possible 
to estimate the total 2 ; ^ 2.5 mass encompassed by IGM 
protocluster candidates. While a full inversion of the 
underlying matter density field from IGM tomographic 
maps is a challeng ing problem beyond the scope of this 
work (although see Nnsse r & Haehnelt 19991 Pichon et al. 
20011 [Kitaura et al.|2U12p, we can use our simulations for 
the more limited goal of estimating the z = 2.5 mass en¬ 
closed within the protocln ster regions. This exploits the 
tight relationship found by ~ 


Lee et al. 


___ ( |2014a| ) between the 

map flux after smoothing by cr = Mpc Gaussian, 

6^^^ and (5dm4, the underlying DM overdensity smoothed 
to the same scale. Figure shows this relationship from 
our simulation (binned in {lh~^ Mpc)^ voxels), which 
can be fitted with a 2nd-order polynomial: 


W ~ 5.502 {6^^y - 3.681 + 0.947. 


( 4 ) 


We then apply this mapping to the values within 
2 : = 2.5 tomographic map voxels encompassed by the 
< —3.5 protoclnster contour, and then sum 
these DM overdensities over the protoclnster candidate 
(‘pc’) region to define a raw ‘tomographic mass’: 


M 


tomo,raw 



X 8.1 X M©, 


( 5 ) 


where M = 8.1 x M© is the cosmic mean den¬ 

sity of DM and baryons multiplied by the (1 h~^ Mpc)^ 
comoving volume of our map voxels. Figure shows the 
estimated 2 : = 2.5 raw tomographic masses for the pro¬ 
tocluster candidates in our simulation against the ‘true’ 
unsmoothed DM mass within each protoclnster region. 
There is a reasonably tight relationship between Mdm 
and Mtomo,raw 5 but the latter systematically underes¬ 
timates by ^ 20 — 30% the true DM mass over the 
^ 1013-14 range typical of our protoclnster re- 
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Figure 8. Scatter-plot of the smoothed Lya forest tomographic 
map flux in our simulation against underlying DM overdensity 
smoothed on the same 4 h~^ Mpc scale, evaluated in (1 h~^ Mpc)^ 
map voxels. The contours denote the 10th, 20th, 30th, 50th and 
80th percentiles of the distribution, while the red circles denote 
the means within bins of AS^ = 0.02 along with the standard 
deviations as error bars. The red-dotted line shows the 2nd-order 
polynomial fitted to the distribution (Eq. |^. 



Figure 10. Descendant cluster mass at 2 : = 0, as a function of the 
corrected tomographic mass (Eq. in 2 ; = 2.5 protoclusters de¬ 
tected in simulated tomographic maps. Red points show the mean 
and standard deviation in logarithmic bins of 0.4 dex in Mtomo, 
while the dotted red line shows the best-fit power-law through the 
points (Eq. [^. The blue dashed line indicates the val ue co rre- 
sponding to our observed 2 ; 2.45 IGM overdensity (Sec. [ 53 . 



Figure 9. Tomographic mass of protocluster candidates in our 
simulation, compared with the true DM mass within each proto¬ 
cluster region at the same redshift of 2 ; = 2.5. The solid black 
line shows the 1:1 linear slope, while the dotted red line shows the 
best-fit power-law fit to correct the tomographic masses (Eq. [^. 


gions, because it was calibrated to the smoothed DM field 
and not the true DM field (binned in {lh~^ Mpc)^ cells 
in our simulation). Since the protocluster boundaries 
(defined by < —3.5) are still significantly over- 

dense, the smoothing smears DM out of the protocluster 
regions hence causing an underestimate of protocluster 
mass. To get an accurate estimate of the 2 ; = 2.5 DM 
mass enclosed within any protocluster candidates iden¬ 
tified in a tomographic map, we would therefore need 
to correct the raw tomographic masses by the power-law 
found to fit Figure 


Mt 


tomo 


10 


14.21 


M 


( 1014 / 1-1 M© 


tomo,raw 


0.946 


/i-i M©. (6) 


This is the formula we use when we henceforth refer to 
‘tomographic masses’. Since this calibration is based 


on tomographic maps with realistic sightline sampling 
and noise, the scatter in Figure provides a measure 
of the uncertainty in the protocluster mass estimate: 
for example, the uncertainty in the mass estimate for a 
Mtomo(^ = 2.5) = 10^^ h~^ Mq protocluster is atomo ~ 
5 X 10^3 M©. 

As suggested by S15, we also checked if the ^ = 2.5 
tomographic protocluster mass could be used to predict 
the 2 : = 0 descendant mass. This is shown for our simu¬ 
lated protocluster candidates in Figure where we see 
a trend for increasing M{z = 0) with Mtomo? which can 
be fitted with the following power-law: 

Part of the ^ 0.2 dex scatter in this relationship is due 
to the reconstruction error in our simulated tomographic 
maps (since we have realistic pixel noise and sightline 
sampling), but S15 have also shown that the dominant 
contributor to this scatter is in fact intrinsic, i.e. due to 
diversity in the morphologies of cluster progenitors at 
fixed M{z = 0). 

Another quantity shown by S15 to correlate with 2 : = 0 
cluster mass is the IGM absorption peak associated with 
the protocluster candidates in the smoothed tomographic 
maps. Figure pd1 shows this for our simulated protoclus¬ 
ter candidates^Again, we see a relationship between the 
absorption depth ((^|?^/crsm)min and the final z = 0 clus¬ 
ter mass, although there is more scatter than the anal¬ 
ogous plot in S15 (Figure 5 in their paper) since we are 
working with realistic reconstructed tomographic maps 
rather than the ‘true’ absorption field in the simulation. 
We fit the following function to the distribution: 

M{z = 0) ^ h~^ M©. (8) 

This absorption depth provides an alternative to Mtomo 
as an estimator for M{z = 0), although both quantities 
are not completely independent. 
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Figure 11. Relationship between the smoothed absorption peak 
associated with all 2 ; = 2.5 protoclusters in our simulated to¬ 
mographic map, and the final z = 0 cluster mass. The red 
points with error bars show the average M(z = 0) in bins of 
A((5|^/crsm)min = 0.5 and the associated standard deviations, 
while the dotted-blue line is the relationship fitted to the distri¬ 
bution (Eq. [^. The shaded gray region is our S^/asm < —3.5 
criterion usea for protocluster selection in the tomographic maps. 
Blue dashed lines indicate the value corresp ondi ng to the observed 
2 ; 2.45 IGM overdensity in our map (Sec. [!3- 

5 . MAP OVERDENSITIES 

In this section, we analyze our reconstructed Lya forest 
map (Section]^ and look for IGM overdensities, aided by 
what we learned from the analysis of simulations in the 
previous section. We will then also discuss two overden¬ 
sities within our volume identified from external galaxy 
datasets. 


5.1. IGM Overdensities 
5.1.1. z = 2.45 Protocluster 


We now apply the S15 protocluster criterion on our 
COSMOS tomographic map. After smoothing by a cr = 
Ah~^ Mpc Gaussian kernel, we find one large IGM over¬ 
density within our map that clearly satisfies the 6^^ < 
—3.5crsm protocluster threshold (Figured. The absorp¬ 
tion peak is at a redshift of 2 ; = 2.450 ano^approximately 
centered on [RA, Dec]?^ [150.06, 2.31] deg, although there 
is a secondary peak at 2 ; = 2.435 within contour de¬ 
fined by the 6^^ < —3.5crsm threshold; these two peaks 
have smoothed absorption values of 6^^ = [—3.9, —4.2] at 
2 : = [2.435,2.450], respectively, which are also marked in 
Figure 1^ This overdensity is highly extended in both the 
line-of-sight and transverse dimensions: it can be seen in 
Figure]^ to continuously extend from 2 : 2.435 — 2.450 

(Slices #3-5), while Figure [T^ shows line-of-sight map 
projections centered on the two absorption peaks of this 
overdensity {z = 2.435 and 2 : = 2.450), showing that it 
spans up to ^ 10 h~^ Mpc 4pMpc) in the transverse 
direction — the overall morphology of the protocluster 
in the unsmoothed map can be seen in Figure The 
total comoving volume of this protocluster candidate (de¬ 
fined as that within the < —3.5crsm contour) is R ^ 
340h“^Mpc^ {7 h~^ Mpc)^. Assuming pure Hubble 

flow, we evaluate the pairwise separations between all 
voxels within the overdensity to find a maximum lin¬ 
ear extent of Tmax ~ 19 h~^ Mpc, which subtends over 


twice our effective spatial resolution element (defined as 
2.3X our effective reconstruction length of 3.5 h~^ Mpc). 
Applying the ‘tomographic mass’ estimation calibrated 
with our simulations (Equations ^ andj^, the total 
DM mass enclosed within the < — To contour is 

Mdm = A/tomo = (9±4) X 10^^ h~^ Mq, where the uncer¬ 
tainty is estimated from the scatter seen in the simulated 
protoclusters (Figure . 

This IGM overdensity corresponds to a known object: 
D15 recently reported a z = 2.450 protocluster com¬ 
prising of 11 LBGs, several of which we plot as filled 
squares in Figure This protocluster was corrobo¬ 
rated by a more recent spectroscopic study of Lya emit¬ 
ters (LAEs) by G15, which are plotted as filled triangles 
in Eigure [l^ The position of these LBGs and LAEs at 
2 : 2.435 — 2.450 appear spatially well-correlated with 

our Lya over density at the same redshifts. The tomogra¬ 
phy also shows that the elongated transverse configura¬ 
tion of the ^ = 2.450 protocluster galaxies is due to two 
sub-structures at that redshift (Eigure 12 r) connected 
by a lower-density ‘bridge’ well-sampleoDy our sight¬ 
lines. Note that due to edge effects, it is possible that 
the protocluster might extend beyond our map boundary. 
Eigures and [T^ suggest that there could be a bridging 
filament between the 2 : = 2.435 and 2 : = 2.450 lobes of 


the overdensity, although deeper observations with bet¬ 
ter sightline sampling and better spectral S/N would be 
required to verify this. 

We now make an estimate of the descendant (z = 0) 
mass of this protocluster based on its IGM tomography 
signature at 2 : = 2.45. The ^ = 2.45 tomographic mass 
enclosed within its < —3.5 boundaries implies, 

through Equation!^ that it will eventually collapse into 
a M(z = 0) = (3.1 ± 2.5) x 10^^ Mq object, i.e. an 
intermediate-mass Virgo-like cluster. The uncertainty is, 
again, based on the scatter of in the simulated proto¬ 
clusters (Eigure pr] with which we calibrated Equation!^ 
The maximum smoothed absorption associated with the 
protocluster also provides an alternative estimate of the 
2 : ^ 0 masses via Equation Applied to the absorption 
peak of ((^|?^/crsm)min ~ —4.2 seen at 2 : = 2.450, this im¬ 
plies 2 : ^ 0 masses of M ^ (3.4 ± 2.1) x 10^^ Mq. 
While these two estimates are not fully independent since 
Mtorno is calibrated through the map pixel values, it is 
reassuring that they give similar values. 


5.1.2. One or two z = 0 clusters? 

Erom the unsmoothed tomographic map (Eigure 0 or 
T^, the 2 : = 2.45 IGM overdensity is comprised ofTwo 
separate lobes at 2 : ~ 2.435 and ^ ~ 2.450. After 
imposing the smoothed < —3.5 cut for pro¬ 

tocluster selection, the resulting contour is continuous 
but still has two absorption peaks: the more signif¬ 
icant one {Spi/asm = —4.2) at 2 : = 2.450 with the 
other at the 2 : = 2.435 lobe with = —3.9, 

with a comoving separation of ~ 16 h~^ Mpc assum¬ 
ing only Hubble flow in the radial direction. Interest¬ 
ingly, if we split our evaluation of the underlying to¬ 
mographic mass to either side of 2 : = 2.442 (the ap¬ 
proximate saddle point between the two lobes) we find 
that the the 2 : = 2.435 part is actually slightly more 
massive at Mtomo ~ 5 x 10^^ h~^ Mq compared with 
Mtomo ~ 4 X 10 ^^ h~^ Mq at 2 : = 2.450. While this 
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Figure 12. Projections of our tomographic map along the line-of-sight, centered on two lobes of the large IGM overdensity identified 
through the < — 3.5(7sm protocluster-finding threshold advocated by S15. The projected redshift range and line-of-sight thickness is 
labeled above each panel, and were chosen to incorporate a reasonable number of coeval galaxies within the projection. Open circles mark 
the sightlines sampling each projected slice, w hile filled symbol s indicate coeval galaxies within the map: MOSDEF galaxies (inverted 
triangles, K15); zCOSMOS-Deep LBGs (stars, |Lilly et al.|(2007|); LAEs from C15 (triangles); and LBGs from D15 (squares). The black 
contour denotes the region that satisfies our protocluster criterion of S^^/asm < —3.5 after smoothing. 





Figure 13. Still frame from the video in Fig. centered on 
the IGM overdensity a,t z ^ 2.44, with the line-of^ight distance 
(z-axis) increasing towards the right. Arrows indicate the two ab¬ 
sorption peaks {z = 2.435 and 2 : = 2.450) associated with the 
z ^ 2.45 protoclus ter, as well as the 2 ; = 2.48 pro toclu ster from 
[Casey et al.| ( [20T^ at the edge of our map (Section |5.2| ). 

difference is insignificant compared to our uncertainties, 
it does point to the overdensity being comprised of two 
roughly equal-mass portions at ^ 15 h~^ Mpc separa¬ 
tion. One might then wonder whether this protocluster 
candidate will in fact collapse into a single z 0 cluster, 
or two separate ones. 

To provide some insight, we again turn to the sim¬ 
ulated protoclusters detected in tomography. S15 de¬ 


scribed such situations in the protocluster search process, 
in which a single 2 ; = 2.5 IGM protocluster candidate 
is in fact associated with progenitors from two separate 
2 ; = 0 cluster^^ Out of the 81 protoclusters identi¬ 
fied through our simulated (256 Mpc)^ tomographic 
map, we find that 21 25%) could in fact be associated 

with a secondary z = 0 cluster. We noticed that many 
of these ‘binary’ protoclusters appeared elongated (e.g. 
middle panel in Figure so we compute on 

their < —3.5 ooundaries as a rough measure 

of their elongation. This is shown in Figure which 
compares with M{z = 0) for both unitary 

and binary protoclusters. In the lower panels of Fig¬ 
ure we see that the binary fraction of the IGM pro¬ 
tocluster candidates is low 10%) for roughly spher¬ 
ical ^ 1) overdensities but increases with 

and approaches unity at Tmax/h^^^^ ^ 2.5, 
i.e. protocluster candidates that appear highly elongated 
in the smoothed IGM tomographic maps are likely to 
evolve into two separate clusters by 2 ; ^ 0. 

Our 2 ; = 2.45 IGM overdensity has a maximum ex¬ 
tent of I/max ~ 19 h~^ Mpc and an enclosed volume of 
V « 340h“^Mpc^, which leads to ^ 2.7. 

This extreme value is only sampled by a handful of sim¬ 
ulated protocluster candidates in Figure but we see 
that ^ 80% of simulated protoclusters with > 2 

collapsed into two separate 2 : = 0 clusters. This indicates 

These were referred to as “bad-merge” objects by S15, al¬ 
though they correspond to distinct M > 10^^ halos at 2 ; = 0 
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L.ax/V^ 


Figure 14. The elongation of IGM protocluster candidates iden¬ 
tified in the simulated tomographic maps. In the top panel, we 
show the descendant mass M{z = 0) associated with each proto¬ 
cluster candidate as a function of the elongation ratio 
In the cases where there is a second cluster progenitor associated 
with the IGM candidate, the secondary protocluster is indicated 
with a blue circle and connected to the primary with a blue line. 
The middle panel shows the distribution of with the 

subset of binaries in blue, while the bottom panel shows the binary 
fraction at each Lmax/I^^'^^ bin. Vertical red dashed line indicates 
the value from the 2 ; = 2.45 protocluster in our observed map. 

a high probability that the IGM overdensity in our map 
will eventually collapse into two separate z = 0 clusters 
with M{z = 0) ^ 10^^ h~^ Mq, centered on the two ab¬ 
sorption peaks at z = 2.435 and z = 2.450. 

5.1.3. Less Significant IGM Overdensities 

Apart from the significant ^ = 2.450 IGM overdensity 
in our map that fulfils the protocluster criterion defined 
by S15, there are also several smaller ‘hot-^ots’ that can 
be seen in the unsmoothed map (Figures]^ and [^, none 
of which fulfill the SfS < —3.5 asm threshed due to their 
limited size. In Slices #4-6, there are two smaller over¬ 
densities at z ^ 2.23 and z ^ 2.25 with transverse extents 
of < 5h~^ Mpe in the lower part of the map (^perp ~ 
6h~^ Mpe and ^perp ~ 3h~^ Mpe, respectively). These 
over densities have map values of « [—2.1, —1.8] asm 
after smoothing and thus do not satisfy our protoclus¬ 
ter criterion, although there are several galaxies associ¬ 
ated with a weaker extension of the z ~ 2.23 structure 
in Slices #5 and #6. There are another two apparent 
overdensities in Slice #1 at z ~ 2.27 and ^ ~ 2.48 at 
^perp ~ 3h~^ Mpe and ^perp ~ 9h~^ Mpe, respectively. 
Neither of these overdensities {dfS ^ [—1.4,—1.8] asm) 
satisfy our protocluster criterion, although this part of 
the map likely suffers from edge effects. A better charac¬ 
terization of these possible overdensities would require an 
extention of the map area beyond its current boundaries. 

5.2. Galaxy Overdensities 

Apart from the overdensity in the Lya absorption map, 
we now discuss two galaxy overdensities from external 


data sets. 


5.2.1. z = 2.30 

During our comparison of the map overdensities and 
coeval galaxies, we noticed a compact overdensity of 
MOSDEF galaxies at 2 ; ^ 2.30 centered on [RA, Dec] 
= [150.06, 2.20]deg (pink circles in Slices #3-4 of Fig¬ 
ure]^; the line-of-sight projection centered on this red- 
shift IS also shown in[T^). This is comprised of 6 galax¬ 
ies within ^ 2 h~^ Mpe (or 1 pMpc) on the sky, and a 
line-of-sight extent of Az ~ 0.006 or Av ~ 560 km s~^, 
corresponding to ~ 6h~^ Mpe (~ 2.6pMpc) assuming 
Hubble flow. If these galaxies fairly trace the underlying 
population, then within the subtended comoving volume 
the number density is pgai ^ 0.3h^Mpc“^, which trans¬ 
lates to an overdensity of Pgai/(/>gai) ^ 40 relative to a 
mean number density of (pgai) ^ 7 x 10“^ Mpc“^, es- 
timated from 2 < z < 3 sources with H < 24.5 in the 
Santini et ahj (|2Q15 ) GANDELS/GOOD-S catalog (see 
Lee et al.|2013b| for a similarly compact z 3.7 overden- 
sity). Within this structure, the quartet of galaxies at 
[RA, Dec] ~ [150.058, 2.200]deg is even more compact: 
with a sky footprint of ^ 0.17arcmin^ and Az ~ 0.003, 
this implies Pgai/(Ppi) ^ 100. 

In Eignres |4| and we do not see any large absorp¬ 
tion decrement associated with these galaxies, although 
the sightline sampling could be better. These galaxies 
occupy map values of 6^^ ^ —0.1 or ^ —0.02 after 
smoothing, therefore based on the calibration between 
smoothed flux and DM overdensity (Equation]^, this 
suggests, at face value, that the MOSDEE galaxies oc¬ 
cupy a region with p/p ^ 1, i.e. the mean density of the 
Universe. This is consistent with the ZEOURGE distri¬ 
bution of K < 2 4.8 galaxies with ph otometric redshifts 
2.1 < 2 : < 2.3 (jSpitler et al.||2012|) in the same field, 
which appears to show roughly mean density at the po¬ 
sition of these MOSDEE galaxies, although 2 : = 2.300 is 
at the edge of the redshift range in the ZEOURGE map. 
Note that since their map is smoothed over a much larger 
radial window {Az = 0.2) than the Az ~ 0.006 line-of- 
sight extent of the MOSDEE overdensity, the ZEOURGE 
map is not in conflict with the presence of the MOSDEE 
overdensity. 

In Eigure we indicate the SfL/ag^ value associ¬ 
ated with this overdensity after smoothing with a a = 
Mpe Gaussian, i.e. the quantity we used to search 
for protoclusters. This corresponds to (^|#/asm ~ —0.4, 
which is far from our (5|#/asm < —3.5 protocluster 
threshold and at a value characteristic of typical regions 
of the map. In comparison with the simulations we see 
that none of the 425 protoclusters within the volume have 
such a weak absorption despite very diverse samplings 
of sightline and spectral S/N across these protoclusters. 
In other words, even the worst-sampled protoclusters 
within our simulated tomographic reconstructions have 
stronger large-scale IGM absorption than this overden¬ 
sity of galaxies. Eurthermore, as can be seen in the bot¬ 
tom panel of Eigure even protoclusters that do not 
satisfy our < —3.5 asm IGM detection threshold still 
exhibit moderate absorption over large ( r\j 10 h ^ Mpe) 
scales which we do not see in the vicinity of these MOS¬ 
DEE galaxies. 

One possibility is that galaxy formation feedback has 
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IGM Tomography of 2 ; = 2.45 Protocluster in COSMOS 
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Figure 15. Same as Fig. |12| but focused on two galaxy overdensities within our map. Panel (a) is focused on a compact overde nsity of 
MOSDEF g alaxies (upside-down triangles), while (b) shows our partial coverage of the 2 : 2.47 protocluster recently reported by |Casey| 

|et al.|(|2015|. Open circles indicate our sightline samplin g, while filled sym bols indicate coeval gala xies within the map : MOSDEF galaxies 
\mverted triangles, K15); zCOSMOS-Deep LBGs (stars, jLilly et al.|20Q7| >, and a DSFG (diamond, |Gasey et al.|2Ql^ . 


suppressed IGM opacity in the immediate vicinity of 
these galaxies on scales of our r\j 3.5 h ^ Mpc effec¬ 
tive smoothing of our tomographic re construction. This 
could a similar effect observed by [Adelberger et ah 


(2005), who claimed a detection of elevated Lya trans- 
mission within ^ 0.bh~^ Mpc of ^ ^ 3 LBGs due to 
galactic winds, althoug h this effect was not seen in hydro- 
dyna mical simulations ( Kollmeier et al.||2006 Viel et al 


2013) 


2UT0 


, nor did subsequen t observations (jbteidel et al 
(Crighton et al.||2011[ ) manage to reproduce this re- 
Hov ' ^ 


suit. However, due to the compactness (mean 3D comov¬ 
ing separation ^ 0.6 h~^ Mpc) of this z = 2.30 galaxy 
overdensity it is possible that the effect of feedback 
has been enhanced beyond that observable in individ¬ 
ual LBGs. Another possibility is that a hot (T > 10^ K) 
intra-cluster medium (IGM) has already formed from the 
virialization of such a compact overdensity, si milar to 


that observ ed in X-rays for a 2 ; = 2.07 cluster (Gobat 
|et al.|[2MT] ). The presence of a hot IGM might also sup¬ 
press the Lya forest absorption on scales of ^ lh~^ Mpc 
(assuming a virial mass of M ^ 10^^-^ h~^ Mq). 

It is possible that a better sampling could reveal a 
stronger overdensity in the immediate vicinity (within 

< 1 — 2Mpc) of this galaxy overdensity, but there is lit¬ 

tle strong absorption in the neighboring sightlines we do 
have, which constrains the transverse extent of the as¬ 
sociated IGM absorption to < Mpc comoving, or 

< 2 Mpc physical. This lack of large-scale overdensity 
suggests that this association of galaxies is unlikely to 
grow into a massive galaxy cluster by ^ ^ 0 (Chiang et al 


plying a very conservative toy galactic feedback model 
on our simulated tomographic maps: for all map voxels 
within Mpc (approximately the virial radius of 

^ 2 X 10^^ h~^ Mq halo) of the absorption peak associ¬ 
ated with each protocluster, we set F = 1 or Sf = 0.25 
(since we adopted {F) = 0.8), i.e. zero Lya absorption, 
and then smooth the maps by a cr = Ah~^ Mpc Gaus¬ 
sian kernel and evaluate the minimum within 

Mpc of each protocluster center as before. The 
resulting distribution is shown as the red dot-dashed 
histogram in Figure which is shifted towards higher 
compared to the original protocluster distribu¬ 
tion as expected, but only by a small amount insufficient 
to match the low map absorption associated with the 
2 : = 2.300 MOSDEF overdensity. 

While it is possible that we have been unlucky with our 
sightline sampling and there is in fact an extended IGM 
overdensity that just happened to avoid our sightlines, 
the current data implies that regardless of any possible 
effects from galactic feedback or a hot IGM on scales of < 
I h~^ Mpc, this 2 : = 2.300 overdensity of galaxies appears 
unlikely to evolve into aM>I0^^h“^ Mq galaxy cluster 
by Z 0 , and could instead b e a precursor of a compact 


similarly com pact 2 : = T.8 LBG overdensity found by|Lee| 


galaxy group (Hickson 


iMTfc. This contrasts with the 


et al. (2013b), in which follow -up observations revealec 
a large-scale LAE overdensity ( Lee et aL]|20I4c ). 


5.2.2. z = 2.48 


Recently, 


2013) due to the lack of available material to accrete from that was initially identihed as a close association of dusty 


Gasey et ahl (l20I5|) reported a protocluster 


its large-scale environment. We demonstrate this by ap- star-forming galaxies (DSEGs), but was also shown to 
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correspond to an LEG overdensity. While the reported 
center of this overdensity lies outside our map volume, 
several of the protocluster members fall just within the 
edge of our map boundary. Their projected transverse 
positions are juxtaposed with our map in Figure [T^. 
Encouragingly, there is an IGM overdensity associated 
with these galaxies right at the edge of our map, which 
is sampled by two sightlines: one in the top-left corner 
of our map area (Figure |T^) and another just outside 
the map boundary (sighthne 15175 in Figure]^. Glearly, 
more data would be required to carry out a aetailed in¬ 
vestigation of this structure. 


large IGM overdensity associated with these galax¬ 
ies; rather, they occupy a region of approximately 
mean absorption on scales of several Mpc. While it 
is possible that galaxy feedback or a hot IGM has 
suppressed the Lya absorption in the immediate 
vicinity (< lh~^ Mpc) of this overdensity, the lack 
of an extended r > Mpc overdensity implies 

that they are unlikely to grow into a massive cluster 
by Z 0. However, the current sighthne configu¬ 
ration does not sample this region well, and more 
observations will be needed to better characterize 
the environment of this overdensity. 


6. GONGLUSION 

In this paper, we describe Lya forest tomographic re¬ 
constructions of the 2.2 < 2 : < 2.5 Lya forest from 58 
background LEG and QSO spectra within a ^ 12' x 14' 
area of the GOSMOS field, which maps the IGM large- 
scale structure with an effective smoothing scale of ^ 
3.5 Mpc over a comoving volume of F ~ 5.8 x 
10 ^ 

Our main findings can be summarized as follows: 

• We compared our map with 61 coeval galaxies 
with known spectroscopic redshifts from other sur¬ 
veys. These galaxies preferentially occupy highly- 
absorbed regions of the map at high statistical sig¬ 
nificance; this skew is even stronger with the sub¬ 
sample of 31 galaxy redshifts measured with NIR 
nebular lines from MOSDEF (K15), which suffer 
less line-of-sight positional error that could scatter 
galaxies into underdense map regions. Future to¬ 
mographic maps encompassing hundreds of coeval 
galaxies will allow an investigation of 2 : ^ 2 — 3 
galaxies as a function of their large-scale IGM en¬ 
vironment. 


• After applying the smoothed flux threshold advo¬ 
cated by SI5 for finding galaxy protoclusters on 
Lya forest tomographic maps, we find one signif¬ 
icant and extended (Tmax ~ 19 Mpc) IGM 
overdensity at 2 : « 2.44 that is as sociated with 
a gala xy pr otocluster repo r ted by [Diener et al. 
(2015) and Ghiang et al. ( |2Q15 ). Within the 
^ 340 h~^ Mpc^ volume of this over density, we 
estimate an enclosed mass of = 2.5) = 

(9 ±4) X 10^^ h~^ Mq. In comparison with the dis¬ 
tributions of size and IGM absorption depth of sim¬ 
ulated protoclusters, we argue that this object will 
likely collapse into a M{z = 0) ^ 3 x 10^^ h~^ Mq 
cluster, although from its highly elongated mor- 
phology ^2.7) there is a high proba¬ 

bility that it will collapse into two separate clusters 
by Z 0. Deeper IGM tomographic observations 
to fully characterize the morphology of this over¬ 
density, along with detailed modeling incorporat¬ 
ing the associated member galaxies, could allow us 
better distinguish between these scenarios. 


These observations were from the pilot phase of the 
upcoming GLAMATO survey, which is aimed at map¬ 
ping the 2 : ^ 2 — 3 IGM across a ^ 1 deg^ area of the 
GOSMOS field. This pilot study validates the ability of 
Lya forest tomographic reconstructions to study large- 
scale structure at these unprecedentedly high redshifts, 
particularly protocluster overdensities that will eventu¬ 
ally collapse into massive galaxy clusters by 2 : ^ 0. A 


1 deg^ survey over 2.2 < ^ < 2.5 would yield comov¬ 
ing volume of F ^ 10^ h~^ Mpc would yield ^ 5 pro¬ 
tocluster detections with M{z = 0) > 10^^ h~^ Mq of 
which ^ 3 would be progenitors of Virgo-like clusters 
{M{z = 0) > 3 X 10^^ h~^ Mq or larger; see also S15). 

Apart from identifying protoclusters with blind tomo¬ 
graphic surveys such as GLAMATO, the insights from 
this paper suggest that it would also be profitable to 
carry out similar observations targeted at known or sus¬ 
pected protoclusters at 2 : ^ 2 — 3. Indeed, over such 
limited fields it could be worthwhile to pursue higher- 
S/N integrations on higher area-densities of background 
sources than we have achieved here, which would allow 
reco nstructions on sm aller scales with reduced map noise 
(see|Lee et al.||20I4a| for a detailed discussion). In syn¬ 
ergy witlTTpectroscopic confirmation of coeval member 
galaxies and boosted by detailed modeling with simula¬ 
tions like those used here, such observations could reveal 
the likely assembly mechanism of individual protoclus¬ 
ters, which would allow workers to build up a global pic¬ 
ture of massive cluster formation from high-redshifts to 


2 : ^ 0. 
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